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We study the hadron-quark phase transition in the interior of neutron stars (NS). For 
the hadronic sector, we use a microscopic equation of state (EOS) involving nucleons 
and hyperons derived within the Brueckner-Hartree-Fock approach. For the quark sector, 
we employ the MIT bag model, as well as the Nambu-Jona-Lasinio (NJL) and the Color 
Dielectric (CD) models, and find that the NS maximum masses lie in the interval between 
1.5 and 1.8 solar masses. 

1. Introduction 

The appearence of quark matter in the interior of massive neutron stars is one of the 
main issues in the physics of these compact objects [ Ij. Calculations of NS structure, 
based on a microscopic nucleonic equation of state, indicate that for the heaviest NS, 
close to the maximum mass (about two solar masses), the central particle density reaches 
values larger than 1/fm 3 . In this density range, it can be expected that the nucleons start 
to loose their identity, and quark degrees of freedom are excited at a macroscopic level. 

The value of the maximum mass of NS is probably one of the physical quantities 
that are most sensitive to the presence of quark matter in the core. Unfortunately, the 
quark matter EOS is poorly known at zero temperature and at the high baryonic density 
appropriate for NS. One has, therefore, to rely on models of quark matter, which contain a 
high degree of uncertainty. In this paper we use a definite nucleonic EOS, which has been 
developed on the basis of the Brueckner many-body theory, and three different models 
for the quark EOS, respectively the MIT bag model, the Nambu-Jona-Lasinio and the 
Color Dielectric models. We compare the predictions of different models, and estimate 
the uncertainty of the results for the NS structure and mass. 

2. EOS of nuclear matter 

Over the last two decades the increasing interest for the equation of state (EOS) of 
nuclear matter has stimulated a great deal of theoretical activity. Phenomenological and 
microscopic models of the EOS have been developed along parallel lines with complemen- 
tary roles. The latter ones include nonrelativistic Brueckner-Hartree-Fock (BHF) theory 
[12] and its relativistic counterpart, the Dirac-Brueckner (DB) theory [HJj, the nonrel- 
ativistic variational approach [ Hj, and more recently the chiral perturbation theory [ 
Ej. In these approaches the parameters of the interaction are fixed by the experimental 
nucleon-nucleon and/or nucleon-meson scattering data. We have calculated the nucleonic 
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equation of state of nuclear matter within the BHF theory. As in all non-relativistic 
many-body approaches based only on two-body forces, the EOS derived in the BHF the- 
ory fails to reproduce some nuclear properties, such as the binding energy of light nuclei, 
and the saturation point of nuclear matter. The usual way of correcting this drawback is 
the inclusion of three-body forces (TBF). In the framework of the Brueckner theory, we 
have adopted two classes of TBF, i.e. a microscopic force [E], based on meson-exchange 
mechanisms, and the phenomeno logical Urbana model [Cj, widely used in variational cal- 
culations of finite nuclei and nuclear matter For details, the reader is referred to Ref.[ 
|H]. We have extended the BHF approach in a fully microscopic and self-consistent way, 
in order to describe nuclear matter containing also hyperons [E]. We have found rather 
low hyperon onset densities of about 2 to 3 times normal nuclear matter density for the 
appearance of the S~ and A hyperons. (Other hyperons do not appear in the matter). 

In order to study the neutron star structure, we have to calculate the composition and 
the EOS of cold, catalyzed matter, by requiring that the neutron star contains charge 
neutral matter consisting of neutrons, protons, hyperons, and leptons (e~, fi~) in beta 
equilibrium. Then we compute the composition and the EOS in the standard way [IT]I10|. 
i.e. by solving the equations for beta-equilibrium, charge neutrality and baryon number 
conservation. The inclusion of hyperons produces an EOS which is much softer than the 
purely nucleonic case. As a consequence, the maximum mass for neutron stars turns out 
to be less than 1.3 solar masses [|H|, which is below the observational limit of 1.44 solar 
masses [ITT]. 

3. Quark matter 

The current theoretical description of quark matter is burdened with large uncertainties, 
and for the time being we can only resort to phenomenological models for EOS, and try 
to constrain them as well as possible by the few experimental information on high density 
baryonic matter. One of these constraints is the phenomenological observation that in 
heavy ion collisions at intermediate energies (10 MeV/A < E / 'A < 200 MeV/A) no 
evidence for a transition to a quark-gluon plasma has been found up to about 3 times 
the saturation density p . We have taken this constraint in due consideration, and used 
an extended MIT bag model [ H2] (including the possibility of a density dependent bag 
"constant") and the color dielectric model [EJ], both compatible with this condition [ 
ITT] . For completeness, we have also used the Nambu-Jona-Lasinio model [IT3]. For the 
description of a pure quark phase inside the neutron star, we have solved the equilibrium 
equations for the chemical potentials of the different quark species, i.e. (u, d, s), along with 
the charge neutrality condition, and the total baryon number conservation. Hence, we 
have determined the composition and the pressure of the quark phase. In order to study 
the hadron-quark phase transition in neutron stars, one has to perform the Glendenning 
construction [HB], by imposing that the pressure be the same in the two phases to ensure 
mechanical stability, while the chemical potentials of the different species are related to 
each other satisfying beta stability. This procedure yields an EOS for the pure hadron 
phase, the mixed phase, and the pure quark matter region. We have adopted a simplified 
method, by demanding a sharp phase transition and performing the Maxwell construction. 
We have found that the phase transition in the extended MIT bag model takes place at a 



The hadron-quark phase transition in neutron stars 



2.4 
2.1 
1.8 h 
1.5 
1.2 h 
0.9 
0.6 h 
0.3 - 




i — • — i — < — i — > — i — <— 

CDM 

MIT, B(p) 

NJL \ 

MIT, B=90 X 



\ 



_l 1 I 1 l_ 



/ 

/ 



V 



/ / / / 
// / / 

// / 

// / ' 

// / ' 

/ / / — ; 



zr 



_i i i i i i i i i i i i i i i i ■_ 



8 9 10 11 12 13 1 
R (km) 



40 0.4 0.8 1.2 1.6 

P c (fm~ 3 ) 



Figure 1. The gravitational mass (in units of the solar mass M ) versus the radius (left 
panel) and the central energy density (right panel). See text for details. 



large baryon density, p ~ 0.6 fm -3 , and at larger baryon density in the NJL model [ IT5] . 
On the contrary, the transition density in the CD model is p ~ 0.05 fm~ 3 . This implies 
a large difference in the structure of neutron stars. In fact, whereas stars built with the 
CD model have at most a mixed phase at low density and a pure quark core at higher 
density, the ones obtained with the MIT bag model contain a hadronic phase, followed 
by a mixed phase and a pure quark interior. The scenario is again different within the 
Nambu-Jona-Lasinio model, where at most a mixed phase is present, but no pure quark 
phase. 

4. Neutron star structure 

We assume that a neutron star is a spherically symmetric distribution of mass in hydro- 
static equilibrium. The equilibrium configurations are obtained by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equations [ for the pressure P and the enclosed mass 
m, 

dp Gm (e + p)(l + 47rr 3 p/m) 

dr r 2 1 — 2Gm/r 



(1) 



dm g , , 

- = 4,A, (2) 



with the newly constructed EOS for the charge neutral and beta-stable case as input, 
supplemented by the EOS of the crust [H]. The solutions provide information on the 
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interior structure of a star, as well as the mass-radius relation, M(R). The results are 
shown in Fig. [TJ displaying mass-radius (left panel) and mass-central density relations 
(right panel). The dashed lines represent the calculation for beta-stable quark matter 
with the CDM, whereas the dotted and dot-dashed lines denote the results obtained with 
the MIT bag model (respectively for a constant bag constant B = 90 MeV fm -3 and a 
density dependent one). The long dashed line represents the calculations obtained within 
the NJL model. We observe that the values of the maximum mass depend on the EOS 
chosen for describing quark matter, and lie between 1.5 and 1.97 solar masses. We notice 
that the inclusion of the color superconductivity in the quark matter EOS built with the 
NJL model decreases the value of the maximum mass down to 1.77 M & thus keeping 
the neutron star maximum mass well below two solar masses. Moreover, neutron stars 
built with the CDM and NJL models are characterized by a larger radius and a smaller 
central density, whereas neutron stars constructed with the MIT bag model are more 
compact, since they contain quark matter of higher density. 

In conclusion, the experimental observation of a very heavy (M > 1.8M Q ) neutron 
star would suggest that either serious problems are present for the current theoretical 
modelling of the high-density phase of nuclear matter, or that the assumptions about the 
phase transition between hadron and quark phase are substantially wrong. In both cases, 
one can expect a well defined hint on the high density nuclear matter EOS. 
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